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ABSTRACT 

We present spectra of the DAB white dwarf PG 1115+166. Radial velocity measure- 
ments of the Balmer lines and the Hel 6678A line show that this is a binary white 
dwarf with a period of 30.09d in which the Balmer lines move in anti-phase to the 
Hel line, i.e., PG 1115+166 is a DA+DB binary. The minimum masses of the stars are 
M DA = O.43±O.15M and M DB = 0.52±0.12M Q . The separation of the stars is about 
45Rq, which is much smaller than a typical AGB progenitor of a white dwarf, imply- 
ing that there has been at least one common envelope phase in this binary. Indeed, 
it is possible that this binary may have suffered up to three mass transfer episodes - 
two associated with the red giant phase prior to the formation of each white dwarf 
and a third associated with the "born-again" red giant phase of the DB white dwarf. 
PG 1115+166 has the longest orbital period of any post common envelope white dwarf 

- white dwarf binary found to date. Published models for the formation of white dwarf 

- white dwarf binaries do not predict any white dwarfs with the combination of long 
orbital period and high mass found in PG 1115+166. We conclude that PG 1115+166 
is a key object for testing models of binary star evolution and it may also be a key 
object for our understanding of the formation of DB white dwarfs. We outline the ob- 
servational tests which can be applied to scenarios for the formation of PG 1115+166 
and apply them to the simplest case of a single common envelope phase. This suggests 
that some part of the internal energy stored in the envelope of the AGB star, e.g., as 
ionized hydrogen, may have contributed to the ejection of the common envelope, but 
there are several unanswered questions concerning this simple scenario. 
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1 INTRODUCTION 

The chemical composition of the atmosphere of a white 
dwarf star reflects a balance between gravitational settling, 
which tends to produce a pure hydrogen atmosphere, and 
processes which "pollute" the atmosphere such as radiative 
levitation, convection and accretion (Fontaine & Wesemael 
1987). This balance will change during the lifetime of a white 
dwarf as it cools and this is seen as a change in the rela- 
tive numbers of hydrogen-rich (DA) and helium rich (non- 
DA) white dwarfs as a function of effective temperature. In 
white dwarfs with effective temperatures T c g ~12 - 60 000K, 
gravitational settling dominates so about 80 percent of these 
white dwarfs have pure or nearly-pure hydrogen atmo- 
spheres. White dwarf stars with helium-rich atmospheres 
may appear as DO stars with He II lines (T e ff Jj45 000K) 
or DB stars with Hel lines (30 000K £ T efi ^12 000K). The 
absence of non-DA white dwarfs with effective temperatures 
between 30 000K and 45 000K is known as the DB gap. The 



favoured explanation for the DB gap is that as the DO white 
dwarf cools to 45 000K, small amounts of hydrogen hidden 
in the atmosphere rise to the surface, eventually masking the 
underlying helium. A hydrogen layer mass of just 10" 14 M Q 
is enough to make the white dwarf look like a DA. Later, 
at ~ 28 000K, the convection zone in the underlying helium 
layer reaches the photosphere, mixing the hydrogen and he- 
lium so the white dwarf re-appears as a helium-rich, DB 
white dwarf. Those white dwarfs with thicker hydrogen lay- 
ers remain as DA stars in this temperature range. In this 
scenario, ~ 20 percent of DA white dwarfs within the DB 
gap have thin hydrogen layers which are sufficient to mask 
the underlying helium. 

The reason why ~ 20 percent of white dwarfs are born 
with very thin hydrogen envelopes is thought to be related 
to the born-again red giant phenomenon in which some stars 
form hot white dwarfs with thick hydrogen envelopes which 
undergo a late thermal pulse. This is caused by re-ignition of 
the helium shell which results in the stars rapidly expanding 
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to red giant dimensions, e.g., Sakurai's Object (Duerbeck 
2000). The remaining hydrogen is quickly mixed into the 
helium layers and is burnt (Iben et al. 1983). 

Those few white dwarfs which show a mixture of hydro- 
gen and helium in their atmosphere are particularly interest- 
ing in the light of this debate. About 20 percent of DB stars 
show traces of hydrogen in their atmospheres and are classi- 
fied as DBA (Shipman, Liebert & Green 1987). Much rarer 
are stars which show both strong hydrogen lines and strong 
helium lines. These are classified as DAB stars and only six 
of them have been identified. The prototype, GD 323, ap- 
pears at the cool edge of the DB gap (T e ff = 28 750K) and 
appears to have a very thin hydrogen layer (10 _17 Mq) al- 
though no model can yet explain all the features of this pe- 
culiar object satisfactorily (Koester, Liebert & Saffer 1994). 
MCT 0128-3846 and MOT 0453-2933 appear to be bina- 
ries with a DA and a DB or DBA component (Wesemael 
et al., 1994). Holberg, Kidder & Wesemael (1990) detected 
a weak HeI447lA feature in G104-27 (WD 0612+177) but 
this was not confirmed by subsequent observations (Kidder 
et al. 1992). HS 0209+0832 is a DAB star in the middle of 
the DB gap which shows strong Hel lines whose strength 
varies from year-to-year (Heber et al. 1997) which may be 
explained by accretion onto the white dwarf from an in- 
homogeneous interstellar cloud (Wolff et al. 2000) . Finally, 
PG 1115+166 was first reported as a DAB star by Burleigh 
et al. (2001). They note that the Ha line shows a variable 
radial velocity but were unable to determine the orbital pe- 
riod. Bergeron & Liebert (2002) have analysed the optical 
spectrum of PG 1115+166 and found an excellent fit to the 
spectrum by assuming it is an unresolved DA+DB double 
degenerate binary. Their fit to the observed spectrum used 
a DA star with T cff = 22090K, log g = 8.12 combined with a 
helium-line DB star with T off = 16210K, logg = 8.19. They 
estimate that both stars have masses of about O.7M0. 

Those DAB stars which are binaries may be particu- 
larly interesting as they may be examples of stars which 
have undergone two mass transfer episodes, e.g., a common 
envelope phase (Iben & Livio 1993). The number of such 
close white dwarf pairs has increased rapidly in recent years. 
The majority of these binaries were identified as low mass 
white dwarfs (^ O.45M0) from analysis of their Balmer lines. 
White dwarfs of such low mass are not formed in the stan- 
dard picture of single star evolution within the lifetime of 
our Galaxy. It now appears that the majority of these bina- 
ries are the result of a common envelope phase in a star on 
the first giant branch (Marsh, Dhillon & Duck 1995). The 
properties of the binaries identified to-date, e.g., their orbital 
periods, masses and mass ratios, are strong tests of models 
for the formation of close binary white dwarfs and mod- 
els for interacting binary stars in general (Nelemans et al. 
2000, Nelemans et al. 2001). Identifying DAB stars may be 
a promising way to find more massive close binary white 
dwarf stars and so extend the range of masses over which 
observations can be used to test models of interacting binary 
stars. 

In this paper we report radial velocity measurements of 
the Balmer lines and Hel lines of PG 1115+166. The Balmer 
lines show that this star is a binary with a period of 30.09d 
and the helium lines move in anti-phase to the hydrogen 
lines, so this is star is a DA+DB binary white dwarf. We 
also consider the formation of PG 1115+166. 



2 OBSERVATIONS AND REDUCTIONS 

PG 1115+166 was first identified as a DAB star in Febru- 
ary 1996 as part of a survey to find low mass white dwarfs 
(Moran 1999). It was re-observed several times thereafter as 
part of the programme at Southampton to find and char- 
acterize close binary white dwarfs. The Leicester group ob- 
served PG 1115+166 independently in April 1999 as a part 
of a programme to find white dwarfs with helium-rich at- 
mospheres in the DB gap by observing hot white dwarfs 
undetected by EUV surveys and also identified it as a DAB 
star (Burleigh et al. 2001) 

Observations of PG 1115+166 have been obtained with 
the intermediate dispersion spectrograph (IDS) on the 2.5m 
Isaac Newton Telescope (INT) and the dual-beam ISIS spec- 
trograph on the 4.2m William Herschel Telescope (WHT) 
both on the Island of La Palma. A variety of gratings have 
been used to cover a variety of wavelengths at different 
times, as detailed in Table [j]. The exposure times used were 
typically 10-30 minutes. Spectra of an arc lamp are taken 
before and after each target spectrum with the telescope 
tracking the star. None of the CCDs used showed any struc- 
ture in unexposed images, so a constant bias level deter- 
mined from a clipped-mean value in the over-scan region 
was subtracted from all the images. Sensitivity variations 
were removed using observations of a tungsten calibration 
lamp. 

Extraction of the spectra from the images was per- 
formed automatically using optimal extraction to maximize 
the signal-to- noise of the resulting spectra (Marsh 1989). 
The arcs associated with each stellar spectrum were ex- 
tracted using the profile determined for the stellar image 
to avoid possible systematic errors due to tilted arc lines. 
The wavelength scale was determined from a polynomial fit 
to measured arc line positions and the wavelength of the 
target spectra interpolated from the calibration established 
from the bracketing arc spectra. Uncertainties on every data 
point calculated from photon statistics are rigorously prop- 
agated through every stage of the data reduction. The nor- 
malized spectrum of PG 1115+166 is shown in Fig. [j] in two 
sections. The spectra used to produce this figure were the 
1999 Dec WHT spectrum, the 1999 Apr WHT spectrum and 
the sum of the three WHT spectra from 2002 Jan. No at- 
tempt has been made to remove the telluric features near 
6900 A and 7300 A. 



3 ANALYSIS 

3.1 Radial velocities. 

To measure the radial velocities we used least-squares fitting 
of a model line profile. This model line profile is the summa- 
tion of one or more Gaussian profiles with different widths 
and depths but with a common central position which varies 
between spectra and is convolved with a Gaussian profile of 
an appropriate width to model the resolution of each spec- 
trum. 

One must be careful to allow for the asymmetry of 
the line due to pressure shifts, particularly when dealing 
with the wings of the higher order Balmer lines (Grabowski, 
Halenka & Madej 1987) or the Hel lines (Beauchamp, We- 
semael & Bergeron 1997). To account for these asymmetries 
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Table 1. Summary of the spectrograph/telescope combinations used to obtain spectra of PG 1115+166 for this study. The slit width 
used in each case is approximately larcsec. 
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Figure 1. The spectrum of PG 1115+166. 
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in the Balmer lines, we add a 3rd or 4th order polynomial 
in the model line profile. We also exclude data from the fit- 
ting process which is more than 1000 km s~ from the rest 
wavelength of the lines measured. This reduces both the ef- 
fect of the pressure shift and the possibility of blending with 
other lines. The radial velocity measurement is determined 
primarily by the centre of the line and the fitting process is 
identical for all lines, so any remaining pressure shifts would 



only be seen as an offset of a few kms -1 between the radial 
velocities measured for different Balmer lines. This is small 
compared to the typical uncertainties in the measured radial 
velocities. 

The pressure shifts in the much weaker Hel lines are 
harder to deal with, so we simply report the radial velocities 
measured using a single Gaussian profile and a low-order 
polynomial for these lines and include a warning here that 



© 1999 RAS, MNRAS 000, 



4 P. F. L. Maxted et al. 



there are systematic shifts between the velocities measured 
for different Hel lines. 

We first normalize the spectra using a polynomial fit to 
the continuum either side of the line of interest. We use a 
least-squares fit to one of the spectra to determine the model 
line profile. A least squares fit of this profile to each spec- 
trum in which the position of the line and the polynomial 
coefficients are the only free parameters gives the measured 
radial velocities in Table ^. We tried three different initial 
values for the radial velocity in the least-squares minimiza- 
tion to ensure we found the optimum value of the radial 
velocity and excluded the results for any spectra where the 
optimum value could not be clearly identified. 

The uncertainties on each radial velocity measurement 
are calculated by propagating the uncertainties on every 
data point in the spectra right through the data reduction 
and analysis. There are inevitably systematic errors in our 
radial velocity measurements, particularly since we are using 
spectra obtained with a variety of instruments. Maxted et al. 
(2000) found that these systematic differences are no more 
than ~ 1 kms -1 for the majority of the spectrographs used 
in this study. To allow for this, we include an additional 
uncertainty per radial velocity measurement equivalent to 
1/40 of the resolution, typically about lkm s~ . This has 
been added in quadrature to the uncertainties given in Ta- 
ble ^. The radial velocity measurements of the Balmer lines 
are clearly variable. Of the Hel line measurements, only the 
measurements of the Hel 6678 lines are sufficiently numer- 
ous and accurate to show any sign of variability. There are 
too few accurate radial velocity measurements of the other 
Hel lines for this variability to be apparent. To quantify this 
we calculated their weighted mean radial velocity, which is 
the best estimate of the radial velocity assuming this quan- 
tity is constant. We then calculated the x 2 statistic for this 
"model", i.e., the goodness-of-fit of a constant to the ob- 
served radial velocities. We can then compare the observed 
value of x 2 with the distribution of x 2 f° r the appropriate 
number of degrees of freedom and find the probability of 
obtaining the observed value of x 2 or higher from random 
fluctuations of constant value, p = 0.023. This suggests that 
the Hel 6678 lines are likely to be variable, but the evidence 
from this test is not conclusive. 

3.2 The orbital period 

We used the 57 radial velocities measured from the Balmer 
lines to establish the orbital period. We used least-squares 
to fit sine waves at 500 000 trial frequencies uniformly dis- 
tributed over the range to 10 cycles per day and recorded 
the value of chi-squared for each trial frequency. Low values 
of chi-squared occur for orbital periods P ps 31d with the 
lowest value of chi-squared being 68.3 for P = 30.0873d. The 
next most likely orbital period is near the one cycle-per-day 
alias of this period, P = 1.0316d, which has a chi-squared 
value of 87.3. 

In order to test the robustness of our period determi- 
nation, we randomly selected 50 of the radial velocity mea- 
surements and fitted sine waves at 100 000 trial frequencies 
uniformly distributed over the range to 1.2 cycles per day 
and recorded the value of the period at which the minimum 
value of chi-squared occurs. We repeated this process 1000 
times and found that in 973 trials this reduced dataset pro- 



duced a minimum value of chi-squared at an orbital period 
near 30.09d. Most of the other trials favoured orbital peri- 
ods near 1 day which, as we have seen, gives a poor fit to 
the complete data set. 

As an alternative test, we used a Monte Carlo simula- 
tion to estimate how often a binary with an orbital period 
of 1.0316d will appear to have an orbital period near 30.09d 
given radial velocity meausurements with the same tempo- 
ral sampling and accuracy as those presented here. We used 
a least-sqaures fit to our data to find the optimum period 
near 1.0316d and the corresponding semi-amplitude of a sine 
wave fit by least squares. We generated 1000 synthetic data 
sets with the same period and semi- amplitude as this fit and 
with the same temporal sampling as our data. We added 
noise to the synthetic data using a psuedo-random normal 
deviates multiplied by the uncertainties associated with each 
datum. For each of these synthetic data sets we fitted sine 
waves at 100 000 trial frequencies uniformly distributed over 
the range 0.01 to 1.2 cycles per day and recorded the value 
of the period at which the minimum value of chi-squared oc- 
curs. We also recorded the best values of chi-squared which 
occur within 10 percent of P = 1.0316d and P = 30.09d. 
Of the 1000 trials, 983 trials resulted in a minimum value of 
chi-squared near 1.0316d and only 17 gave a minimum value 
of chi-squared at a period near 30.09d. Furthermore, in all 
the trials for which P ~ 30.09d gave a lower value of chi- 
squared, the difference in the minimum chi-squared value for 
P ~ 1.0316d and P ~ 30.09d was usually only 1-3 and was 
always less than 8. In summary, this test suggests that there 
is a less than 1/1000 chance of a binary with P = 1.0316d 
being incorrectly identified as a binary with P — 30.09d 
from our data provided the least-squares fits of sine waves 
at the two periods result in a difference in chi-squared of 8 
or more. The difference in chi-squared for our actual data is 
19, so this condition is easily satisfied in this case. 

These simple test both suggests that period determina- 
tion is robust and that we can be confident that the orbital 
period is P » 30.09d. 

3.3 Spectroscopic orbit 

If PG 1115+166 is genuinely a binary DA+DB white dwarf 
we expect that the Hel lines should move in anti-phase 
to the Balmer lines. We investigated this scenario using 
a least-squares fit of a sine wave of the form Vi = 7h + 
Kh sin(27r(T; — To)/P to the measured radial velocities, 
Vi, for all the Balmer lines and a sine wave of the form 
Vi — 7Hci + KhcI sin(27r(Ti — To)/P to the 23 measured ra- 
dial velocities for the Hel 6678 line. We used both data sets 
in a simultaneous fit to all these data, an approach which 
has the merit of using all the available data to determine the 
parameters To and P. The results are given in Table ^ and 
are shown in Fig. |j| We have excluded data for the Balmer 
lines with uncertainties larger than 20 km/s and data for the 
Hel 6678 line with uncertainties larger than 40 km/s from 
Fig. H for clarity. The value of Knel derived is about 4 times 
greater than its uncertainty, so we can be fairly sure that 
the Hel 6678 line moves in anti-phase to the Balmer lines as 
expected for a DA+DB binary. We can certainly rule out 
the possibility that PG 1115+166 is a single DAB star since 
that would require the Hel lines and the Balmer lines to 
move in phase, which is clearly not the case. Also given in 
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Figure 2. Radial velocities measured from the Balmer lines (filled 
symbols) and the HeI6678A line (open symbols) of PG 1115+166 
as a function of orbital phase. The sine waves fit by least-squares 
described in the text arc also shown. 
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Table ^ are the results of a similar least-squares fit including 
the data for all the Hel lines. This reduces the uncertainty 
in the value of Kuei , but we prefer the results of fitting only 
the Hel 6678 line until the systematic errors introduced by 
pressure shifts in Hel lines is better understood. In any case, 
the results of the two least-squares fits are nearly identical. 

The mass ratio implied by the parameters in Table |^ is 
q = A/da/A/db = 0.84 ± 0.21. This is consistent with the 
conclusion of Bergeron & Liebert (2002) that the masses of 
the two stars are similar. We also find Mda sin 3 i = 0.43 ± 
0. 15M , A/db sin 3 i = 0.52 ± 0. 12M Q and a sin i = 40 ± 5R , 
where i is inclination of the binary and a is the separation 
of the stars. The masses measured by Bergeron & Liebert 
suggest that the inclination of the binary is approximately 
60° and that the stars are separated by about 45R©. 



4 THE FORMATION OF PG 1115+166 

The basic mechanism for the formation of most white dwarfs 
is fairly well understood. They are, in general, remnants 
of intermediate mass stars (I-8M0) which have evolved 
through the asymptotic giant branch (AGB) phase. Extreme 
mass loss towards the end of this phase exposes the core of 
the AGB star, which cools to form a white dwarf. Stars 
with larger initial masses are expected to form more mas- 
sive white dwarfs and there is some observational evidence 
for this initial-final mass relation (IFMR; Weidemann 2000) . 
White dwarfs with masses of O.7M0 are expected to form 
from stars with initial masses of about 3M0. 

Models of stars on the AGB are uncertain, but the ra- 
dius of a star towards the end of this phase is several hun- 
dred solar radii (Girardi et al. 2000). This is much larger 
than the current separation of the stars so it appears that 
the orbit of PG 1115+166 is now substantially smaller than 
it was prior to the formation of the white dwarfs. Orbital 
shrinkage in long period binaries is explained by the "com- 
mon envelope" scenario. In this scenario, the expanding red 
giant star comes into contact with its Roche lobe and begins 
to transfer mass to its companion star. This mass transfer 
is highly unstable, so a common envelope forms around the 
companion and the core of the red giant. The drag on the 



Figure 3. PG 1115+166 in the mass-period plane compared to 
the models of Iben, Tutukov & Yungelson (1997). The figure is 
taken from Saffer, Livio & Yungelson (1998) and shows the pre- 
dicted number of observable white dwarf- white dwarf binaries 
as a function of orbital period and mass as a gray scale images. 
We have assumed an uncertainty of 0.05Mq in the mass for this 
figure. 
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companion orbiting inside the common envelope leads to ex- 
tensive mass loss and dramatic shrinkage of the orbit (Iben 
& Livio 1993). We expect the orbit of a post-CE binary to 
be circular, and have assumed this to be so in the case of 
PG 1115+166 throughout this paper. Much higher quality 
data than that presented here would be required to verify 
this assumption. 

Common envelope evolution is poorly understood so 
it is useful to consider the limits imposed on the com- 
mon envelope (CE) phase by the observed properties of 
PG 1115+166. Firstly, we note that published models for 
the population of white dwarf binaries do not predict any 
systems with the combination of high mass and long period 
seen in PG 1115+166. This is shown in Fig. ^|for the models 
of Iben, Tutukov and Yungelson (1997), but the same applies 
to the models of Nelemans (2001) and Han (1998). It is pri- 
marily the combination of high mass and long orbital period 
which is not explained by these models. A full exploration 
of all the possible formation scenarios for PG 1115+166 is 
beyond the scope of this paper, not least because there may 
have been up to three mass transfer episodes - two asso- 
ciated with the red giant phase prior to the formation of 
each white dwarf and a third associated with the born-again 
red giant phase of the DB white dwarf. These mass transfer 
episodes may not have been CE phases but may instead have 
involved stable mass transfer, i.e., an Algol-like phase. In 
this section we simply outline the observational constraints 
which can be applied to these scenarios and apply these con- 
straints to the simplest case of a single CE phase. 

Apart from the present day masses of the white dwarfs 
(~ O.7M0) and their current separation (a/ ~ 45R©), we 
also know the cooling ages of the stars are 6 x 10 7 y for the 
DA white dwarf and 2.2 x 10 8 y for the DB white dwarf 
(Bergeron & Liebert 2002). The cooling ages are the time 
elapsed since the AGB phase. If the progenitors of the white 
dwarfs were formed at the same time, this imposes a limit 
on the initial masses of the progenitor stars as follows. If 
we assume an initial mass for the progenitor of the DB 
star, M[, we can use an appropriate stellar model to find 
the lifetime of this star, n, i.e., the time taken since it 
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Figure 4. The relationship between the initial masses of the pro- 
genitors of the DB star (M[) and the DA star (Mj) implied by the 
difference in cooling ages of (1.6 ±0.1) X 10 8 y based on the mod- 
els of Girardi et al., (2000) for Z=0.019 (solid lines) and Z=0.004 
(dashed lines). 




formed for the star to reach the end of AGB phase. We then 
know that the progenitor of the DA star will take a further 
2.2 x 10 8 y-6 x 10 7 y = 1.6 x 10 8 y to reach the AGB phase. We 
can then use grid of stellar models to find the initial mass, 
M2, of a star whose lifetime is t% = t\ + 1.6 x 10 8 y, i.e, the 
mass of the progenitor of the DA star. The relationship be- 
tween these initial masses is shown in Fig. ^ based on the 
models of Girardi et al. (2000) for Z=0.019 (solid lines) and 
Z=0.004 (dashed lines). Although the stars have not evolved 
as single stars, this argument regarding the cooling age dif- 
ference between the stars remains valid because the AGB 
phase is much shorter than the lifetime of the star which 
precedes it, i.e., the stars spend most of their lives evolving 
as single stars. However, we should be aware that the IFMR 
may not apply to the white dwarfs in PG 1115+166 in the 
sense that a star may produce a white dwarf which is less 
massive than expected if it is involved in a CE phase which 
interrupts the growth of the future white dwarf forming in 
its core. Also, this constraint may not be valid if there has 
been significant accretion onto one of the stars. Similarly, 
if the born-again red giant phase of the DB star was short, 
the DB white dwarf may not have been re-heated to the 
temperatures expected for a white dwarf emerging from its 
first AGB phase. In this case the cooling age does not mea- 
sure either the time since the AGB phase, nor the time since 
the born-again read giant phase, but is somewhere between 
these two timescales. 

The observation that one of the stars is a DB white 
dwarf is also a constraint on possible formation scenarios 
for PG 1115+166. Any feasible scenario must allow for the 
DB white dwarf to undergo a born-again red giant phase or 
include an alternative explanation for how this star came 
to be hydrogen deficient. A stronger constraint, perhaps, 
is how this star avoids accreting hydrogen onto it surface 
despite appearing to be older than its companion which has, 
presumably, lost large quantities of hydrogen rich material. 

With this information we consider the simplest case of 
a single CE phase in which the AGB star progenitor of the 
DA star with mass M2 and radius R2 interacts with the DB 
star when it is close to or exceeds its Roche lobe. Note that 



M2 will be less than M'2 if any mass loss occurs prior to 
the interaction. However, this mass loss cannot have been 
very extensive in this scenario because the DB white dwarf 
cannot have gained more than ^ 10~ 14 Mq so we assume 
M2 = M2. If the initial semi-major axis of the stars' orbit is 
en and the current separation is a/, the change in the orbital 
binding energy is 

. „ GM2M0B GMdaMdb 

A-C-orb = 

2oi 2a f 

Some fraction of this energy qqe is used to expel the enve- 
lope of the AGB star whose binding energy is given by 

GM 2 {M 2 - M DA ) 



E, 



Ai?i 



So from the definition of acE we have 



qceA 



/ GMdaMdb _ GM 2 M OB \ = GM 2 (M 2 - M DA ) 



2a 



f 



2ai 



Ri 



If the orbit prior to the interaction has an eccentricity e, 
then the Ri — rLai(l — e) at the time of interaction, where 
tl is the radius of the Roche lobe relative to the separation 
of the stars. This assumes that the interaction occurs at 
periastron, which is reasonable given that the orbital period 
(~ lOy) is much shorter than the evolutionary timescale of 
an AGB star. 

In Fig. ^ we show lines of constant cuceA as a func- 
tion of M 2 and Ri/(1 - e) for M DA = M DB = O.7M and 
af = 45R©. The maximum radius of an AGB star with an 
initial mass of ~ 3Mq is a few hundred solar radii (Bloecker 
1995), so unless the eccentricity of the orbit prior to the CE 
phase was quite extreme (e <; 0.5) and/or the initial mass 
of the AGB star was lower than suggested by the IFMR 
(« 2M©), we see that qceA <; 1.0. This requires a large 
value of A because «ce < 1 in the absence of other sources 
of energy to eject the envelope. This strong limit on A is a 
consequence of the combination of long orbital period and 
high mass in PG 1115+166. Large values of A for AGB stars 
are predicted by Dewi & Tauris (2000) provided the cal- 
culation of the binding energy of the envelope includes the 
internal energy stored in the envelope, e.g., as ionized hy- 
drogen. This implicitly assumes that this potential energy is 
efficiently converted into the kinetic energy of the envelope 
during the common envelope phase. Several issues arising 
from this scenario remain to be explored, e.g., how does 
the DB star avoid gaining more than 10~ 14 Mq of hydrogen 
during the CE phase, why was there no CE phase when the 
progenitor of the DB star went through its AGB phase and 
where does the born-again red giant phase fit into this sce- 
nario? It is beyond the scope of this paper to tackle these 
questions, but they do demonstrate that our observations of 
PG 1115+166 provide a challenging test of models for binary 
star evolution and the formation of DB white dwarfs. 



5 CONCLUSION 

We have used radial velocity measurements of the Balmer 
lines and Hel lines of PG 1115+166 to show that this star 
is a DA+DB binary white dwarf with an orbital period of 
30.09d. Published models for the formation of white dwarf 
- white dwarf binaries do not predict any white dwarfs 
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Figure 5. The efficiency of the common envelope phase, a^E^i 
as a function of the radius, R2 and the mass M2 of the AGB star 
progenitor of the DA white dwarf in PG 1115+166. Solid lines 
correspond to values of a^X as marked. The effect of a change 
of 0.02Mq in the current masses of the white dwarfs is illustrated 
by the dotted lines either side of the solid lines 
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with the characteristics of PG 1115+166. We consider the 
observational constraints on scenarios for the formation of 
PG 1115+166. We apply these to the simplest case of a sin- 
gle common-envelope phase. We find that this requires that 
some part of the internal energy in the envelope of the AGB 
star contributes to the ejection of the common envelope, 
but several questions regarding this scenario remain to be 
answered. We conclude that PG 1115+166 is a key object 
for testing models of binary star evolution. It may also be 
a key object for our understanding of the formation of DB 
white dwarfs. 
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Table 2. Radial velocities measured from the Balmer lines and 


Table 2 — continued 
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Table 3. Results of least-squares sine 


wave fits to the measured 
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